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ABSTRACT 

We investigate whether a baryonic dark galaxy or 'galaxy without stars' could persist indefinitely 
in the local universe, while remaining stable against star formation. To this end, a simple model has 
been constructed to determine the equilibrium distribution and composition of a gaseous protogalactic 
disk. Specifically, we determine the amount of gas that will transit to a Toomre unstable cold phase 
via the H2 cooling channel in the presence of a UV — X-ray cosmic background radiation field. 

All but one of the models are predicted to become unstable to star formation: we find that, in the 
absence of an internal radiation field, the majority of gas will become Toomre unstable in all putative 
dark galaxies with baryonic masses greater than 10^ Mq, and in at least half of those greater than 10^ 
M0. Moreover, we find that all our model objects would be detectable via Hi line emission, even in the 
case that star formation is potentially avoided. These results are consistent with the non-detection of 
isolated extragalactic Hi clouds with no optical counterpart (galaxies without stars) by the Hi Parkes 
All-Sky Survey. 

Additionally, where star formation is predicted to occur, we determine the minimum interstellar 
radiation field required to restore gravothermal stability, which we then relate to a minimum global 
star formation rate. This leads to the prediction of a previously undocumented relation between Hi 
mass and star formation rate that is observed for a wide variety of dwarf galaxies in the Hi mass 
range 10^ — 10^° Mq. The existence of such a relation strongly supports the notion that the well 
observed population of dwarf galaxies represent the minimum rates of self-regulating star formation 
in the universe. 

Subject headings: Galaxies: Formation — Galaxies: Kinematics and Dynamics — Galaxies: ISM — 
Stars: Formation — Atomic Processes — Molecular Processes 



1. INTRODUCTION 

There has been long-standing speculation that opti- 
cally selected galaxy surveys might be missing a popula- 
tipn of very low surface brightness (LSB) objects l)Disnevl 
^76: Disnev & PhillioDS 1983) • It has recently become 
possible to address this question observationally: the Hi 
Parkes All-Sky Survey (HiPASS) has now mapped the 
entire S outhern sky at ^ 21 cm. Since the HiPASS cat- 
alogue ijMever et al.ll200'^ represents an optically unbi- 
ased sample of neutral, atomic hydrogen (Hi) in the local 
[z < 0.04) universe, it is ideal for the detection of Hl-rich, 
baryonic dark galaxies; 'galaxies without stars'. 

The result of the HiPASS search for dark galaxies 
is emphat ic and su rprising. With few poss i ble ex cep- 
tion s fe.o.lKilborn e t al. (2000); Rvd er et alJ l|200lD : see 
also lMinchin et alTpOOS )]. no isolated, extragalactic Hi 
clouds were discovered for which an opti cal counterpart 
could not be identified l|Dovle et alJl2005D ; it seems there 
are no dark galaxies. 

Prior to HiPASS, the discovery of a number of 
'dim' galaxies — "galaxies where luminous matter is 
only a minor component of the total galaxy mass" 
ijCarignan fc FreemanllToS^ — had led to the expecta- 
tion that star formation (SF) might be suppressed com- 
pletely within some galaxies. The first such examples 
were dwarf irregular galaxies with extended Hi envelopes; 
DDO 154 se rves as the archet ype of this class of gas-rich 
LSB dwarfs l)Krumm fc Burs tein 1984). At the other end 
of the mass scale, the serendipitous discovery of Malin 1 
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l|Bothun et alJl987t lhm )ev fc Bothunll989() . a very mas- 
sive, Hl-dominated LSB disk galaxy, provided a n exam- 
ple of the class of "crouching giants" predicted bv lDisnevI 
Jl976). 

A common feature of Malin 1 and DDO 154 is their 
high Hi mass:light ratios: Mm/ Lb ~ 5 [compared to 
Mm/Lfi ^ 0.1—1 for 'normal' Hl-detecte d LSB galaxies 
l|Fisher fc T uUv 1975': 'Wau gh et al.l l2002')]. In addition 
to the observational problem of determining their abun- 
dance — which HiPASS has begun to address — galaxies 
like DDO 154 and Malin 1 pose the following theoreti- 
cal problem: "Why was there so little gas processed into 
stars" in these systems ijCarignan fc Beaulieul ll989)? 

Are these objects exceptional in the way in which they 
are presently forming stars? How far might this 'dim' 
population extend toward complete darkness? And if 
dark galaxies can exist, is there some mechanism to pre- 
vent their detection by HiPASS? 

In this paper, we consider the viability of dark galax- 
ies in the local universe; that is whether a dark galaxy, 
having formed, can remain dark, or whether it will in- 
evitably 'fight up'. A simple model has been developed 
for the distribution and composition of primordial gas in 
a protogalactic disk that is in chemo-dynamical equilib- 
rium, in the presence of a photodestructive and photo- 
heating UV — X-ray cosmic background radiation (GBR) 
field. The specific goal is to determine where gravother- 
mal instability (and hence SF) is initiated through the H2 
cooling channel. Additionally, where H2 cooling cannot 
be checked by the CBR alone, we place a lower bound 
on the global star formation rate (SFR) in these objects 
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by determining the additional interstellar radiation field 
(ISRF) required to restore thermal balance. 

The discussion proceeds as follows. First, an overview 
of the argument is provided in ^ as well as physical mo- 
tivation for our model for a protogalaxy. In SJ3 the link 
between H2 cooling and Toomre instability is discussed, 
and we advance our specific criterion for SF. We provide 
a technical discussion of the processes by which we de- 
termine the gas distribution and composition in 21 ^'^d 
^ respectively. The second calculation, in which we es- 
timate the global SFR, is then described in ^ After a 
discussion of fiducial parameter choices in the results 
of our modelling procedure are given in !JH1 and discussed 
in S]3 A summary of the salient and novel aspects of this 
work can be found in HIOI 

2. OVERVIEW AND 
SUPPORTING CALCULATIONS 

2.1. Overview of the Argument 

Our primary goal is to determine whether a puta- 
tive dark galaxy is gravothermally stable in its long 
time/steady state configuration. Where this is the case, 
it is conceivable that such an object (assuming that it can 
form) might remain dark indefinitely; otherwise, from 
this point, if not well before, the galaxy will inevitably 
'light up'. 

Whereas rotational support stabilises against the col- 
lapse of larger-scale density fluctuations, small-scale fluc- 
tuations can only be stabilised by thermal pressure (see 
Once the galaxy formation process is complete, the 
level of rotational support becomes fixed; stability is then 
principally determined by temperature, T. In particular, 
as will be shown in S]31 most low-mass protogalactic disks 
are stable at T = 10"* K, while essentially all will become 
unstable if they are able to cool to T '--^ 300 K. 

The final temperature of the gas depends on its compo- 
sition. Internal cooling of the gas proceeds by radiative 
de-excitation of collisionally excited atoms/molecules. 
Whereas Hi cooling processes rapidly become inefficient 
(in the sense that the timescale for cooling becomes 
greater than the age of the universe) as collisional exci- 
tation ceases below ~ 10^ K, II2 cooling remains efficient 
above ~ 300 K. The key to initiating gravothermal in- 
stability and SF in low mass disks is therefore efficient 
H2 cooling; conversely, the prevention of SF requires the 
prevention of H2 cooling. 

Within this paradigm, SF is thus a process that must 
be prevented, rather than initiated: a source of dissociat- 
ing radiation is required to prevent the partial conversion 
of Hi to H2. Within the model, the CBR is the agent for 
dissociation, as well as the source of heating that opposes 
H2 cooling. SF is deemed to be preventable wherever the 
rate of photoheating due to the CBR exceeds the rate of 
H2 cooling within the gas; otherwise SF is inevitable. 

We thus reduce the problem of determining the 
in/stability of the protogalactic disk to determining the 
in/efficiency of H2 cooling: we assume that gravother- 
mal stability implies and is implied by thermal balance 
at T ^ IC* K. (We will be able to explicitly check the con- 
verse assumption — that thermal instability at T = 10^ 
K implies cooling to Toomre unstable temperatures — in 
TO. 21 ) Because we are interested primarily in the situa- 
tion of permanent gravothermal stability, we first assume 
that the gaseous disk is in dynamical and chemical equi- 



librium at T = 10^ K; ie. a long-time configuration at 
the end of Hl-dominated thermal evolution. 

This configuration represents the long-time state of a 
galaxy that has managed, for whatever reason, to avoid 
SF, since the gas is Toomre stable so long as T > 10^ K; 
if, in this configuration, the gas cannot initiate H2 cool- 
ing, this model therefore represents the final state of a 
galaxy without stars. Conversely, where H2 cooling is 
inescapable, the result will be thermal, chemical, and dy- 
namical instabilities, leading necessarily to gravothermal 
collapse and SF. 

In other wo rd s. foll owing authors like 
ICorbeHi Sa,]^^ flflfll. IHaiman. Thnu1 T.oebl 
(119961) . lOh fc HaimanI 1120021) . and ISchavd (12001) 7 
we link SF to the existence of a cold phase within 
the gas, and examine under what conditions the 
transition to the cold p hase is initiated. [See, too, 
lElmegreen fc Parravand l)1994() . but also lYoung et alJ 
l)2003j) .] Note that turbulent effects, which can both 
initiate and prevent SF, are neglected; all discussion 
of this point is deferred until i)9.1l We have thus 
constructed a static problem — we make no attempt to 
track the evolution of the gas at any stage. 

Wherever possible, our assumptions have been tailored 
to minimise the rates of both H2 production and H2 cool- 
ing, so as to make SF as difficult as possible. We thus 
place a firm lower limit on the mass of gas that can exist 
in a steady state without undergoing any SF. 

Where we find that SF cannot be prevented by the 
action of the CBR, we perform a second calculation to 
place a lower limit on the expected SFR. This is done 
by introducing a diffuse ISRF as a second source of pho- 
toheating and photodissociation. Specifically, we deter- 
mine the minimum ISRF intensity required to maintain 
thermal balance at T = IC* K, and so restore gravother- 
mal stability. With some simple assumptions (see 
this ISRF can then be related to a minimum global SFR 
within the model galaxy. 

2.2. The Back of the Envelope Calculation 

The CBR can act to check H2 cooling in two ways: by 
direct photodcstruction of H2 and by overwhelming any 
H2 cooling with photoheating. The cloud's only defence 
is self-shielding; ie. absorption of enough of the CBR by 
the outermost gas to shield inner regions from its influ- 
ence. 

The action of the CBR is thus to prevent SF in a shield- 
ing 'skin' that covers both faces of a protogalactic disk, in 
a manner analogous to photodisso ciation regions around 
UV s ources in our own galaxy l)Hollenbach fc Tielentj 
^1^; SF is only possible in the self-shielded centre. The 
thickness of this skin can be characterised by the critical 
surface density of gas, Scrit, below which the CBR is able 
to keep the gas completely photodissociated. 

We obtain a flrst estimate for Ecrit by noting that for 
the CBR to prevent the production of H2 , the total rate 
of photodcstruction in a column (which is just the flux 
of photodestructivc photons) must exceed the total re- 
combination rate; viz.: 

^ Ai. > n Zrf . (1) 

Here, the flux of destructive photons is approximately the 
product of the number of photons at the threshold fre- 
quency, Jy/hv (where is the CBR flux at the threshold 
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frequency , and the frequency interval over which pho- 
todestructive absorption can occur, Aiy. Similarly, the 
recombination rate in the column is approximately the 
product of the recombination rate per unit volume, n, 
and some characteristic column height, Zd- 

With the approximations S « poZd (where po is, for 
now, a characteristic density), and h « aXiXjPg/ fjl^ (a 
is the rate coefficient of the reaction between species i 
and j that replenishes the cloud's H2, Xi is the relative 
abundance of species i, and p ~ 1.22mp is the mean 
molecular mass of particles in the gas), it is possible to 
rewrite inequality as: 

\ hv aXiXj 2TrG J 
In writing this expression, we have also made use of equa- 
tion (|17|l to relate po and Zd, which introduces the gas 
velocity dispersion, a. 

Using equation |(2J), we can now obtain a quantitative 
estimate for Ecrit, assuming T — 10"^ K (which fixes cr 
and the value of a), chemical equilibrium (which sim- 
ilarly fixes Xi and Xj), and a putative CBR spectrum 
J(z^). To estimate Av, we assume that H2 shielding oc- 
curs between the Lyman- Werner and Lyman-a threshold 
energies, 11.2 — 13.6 eV. With these numbers, we calcu- 
late Ecrit ~ 3 — 4 Mq pc~^, depending on the character 
of the CBR spectrum (see ^7.3(1 . 

In reality, the CBR will produce at least two layers 
of skin: one in which it prevents Hi production [ie. an 
ionised layer) and, internal to that, one in which II2 pro- 
duction is prevented (ie. a neutral, atomic layer). Using 
the same argument as above, we find that the surface 
density of the ionised layer is 1 Mq pc~^. 

We therefore expect that the CBR will completely pre- 
vent SF where the surface density of the protogalactic 
disk is less than several M© pc~^. [For comparison, ob- 
servations are consistent with a fixed SF threshold of 
~ 10 Mq pc"2 (Taylor et al. 1994; see also Martin & 
Kennicutt 2001 and discussion in Schaye 2004).] 

2.3. Characteristic Timescales 

We assume chemo-dynamical equilibrium in our mod- 
elling. Thus, while our results do not rely on any partic- 
ular galaxy formation scenario, they do depend on these 
equilibrium states being accessible to each system within 
its present lifetime. Our model does not apply to form- 
ing galaxies, nor does it attempt to follow the behaviour 
of the gas as it transits from warm to cold. 

In particular, our criterion for SF implicitly assumes 
that enough time is available for the gas to cool from 
T ~ lO'' K to T - 300 K after the onset of H2 coohng. 
In the remainder of this section, we present some simple 
calculations in support of this assumption. 

First, we can evaluate the timescale for cooling 
tc = e/e, where e ~ nkT is the thermal energy den- 
sity {n is the number density of the gas; k is Boltz- 
mann's constant) and e = A is the net cooling 
rate per unit volume. Using the rates of cooling 
listed in Appendix ^ and assuming gas composition 
in chemical equilibrium at T = 10^ K (see Figure 
ini, we find tc ~ 1.5 x 10^^ n^^ tn = 20 Myr where 
tn ~ 13.7 Gyr is the Hubble time and n has units of 
cm^"^. (Close to the galactic plane, n is typically in the 
range 10~^ — 10^ cm~'^; see Figure^]) 



Under these conditions, the Hi and H2 cooling rates are 
approximately equal. In contrast, the same composition 
at 8000 K leads to cooling times for Hi and H2 cooling 
of 2.3 tn and 15 Myr, respectively: H2 cooling 
remains efficient well below the point at which Hl-driven 
thermal evolution effectively ceases. 

As we argue in §5.3, in gas near thermal equilibrium, 
the H2 abundance will always approach equilibrium from 
above; in this case the equilibrium value therefore pro- 
vides a lower limit to the actual H2 abundance. For an H2 
abundance ten times greater than the equilibrium value, 
XH2 ^ 10~^, the cooling time at T = lO'' K is reduced 
by a factor of 30 to ~ 0.67n~^ Myr. Among tidal dwarf 
galaxies, the only observed population of presently form- 
ing galaxies , the H2 abundanc e is typically as high as 
'^0.20 l|Braine et al.ll200H) : the cooling time in this 
case is reduced to just ~ 10"'^ Myr. 

Next, we can make a rough estimate for the 
time needed to produce such an overabundance, 
= xii2 /^ii2 ■ Again assuming equilibrium abun- 
dances at 10"* K, and using the rate for the dominant 
H2 production process (Reaction 12 in Appendix 
XH2 ~ 2 X 10"^'' s~^, we find ~ 1-6 n~-^ Myr to 
reach — 10^^. This value is of the same order of mag- 
nitude as if we were to take a semi-empirical rate for H2 
production on dust grains at T < 10^ K (two-body H2 
production becomes inefficient at these temper atures), 
XH2 ~ 10~" s-i llHollenbach fc Tielenslll999j) . 

Thus, as equilibrium at 10'' K is disturbed, we have 
^Ha < tc, and xb_2 increases above the equilibrium value; 
tc is accordingly diminished. At the same time, t^^ in- 
creases as the temperature falls, until tc -C tn^. At 
this point two-body chemical evolution effectively ceases, 
leaving a 'freeze-out' abundance of H2: time- integrated 
studies of shock-heated {ie. T > 10^ K, x^- « 1) gas 
coolin g cata strophically below 10^ K (jShaoiro fc Kand 
Il987t lOh fc HaimanI l2002t) find that asymptotes 
to a 'universal' freeze-out abundance of order 10"'^ at 
T ~ 4000 K. With the provision that this fre eze-o ut H2 
abundance can survive photodissociation (see ^9.2|l from 
this point on, the gas continues to cool until tc > % at 
T < 300 K. 

Once the warm (T ~ 10"* K) gas initiates H2 cool- 
ing, it will therefore transit to a cold (T ~ 300 K) 
state in < 20 Myr. This should be compared to 
the characteristic timescale for the first stage of SF, the 
free-fall time for a collapsing cloud, tf — {Sn/lQGpy^^ 
~ 52 Myr. Where H2 cooling occurs, the cloud 

will therefore cool on a timescale at least a few times 
faster than the collapse timescale; depending on the ac- 
tual H2 abundance, the factor is more likely > 30. 

3. THE MODEL — I 
THE CRITERION FOR STAR FORMATION 

The standard criterion for the gravotherma l stabil- 
itv o f a rotating disk, the Toomre criterion (^Toomr j 
IT961 . can be obtained from the dispersion relation 
for an axisymmetric density perturbation (ie. a sound 
wave) with wavelength £ and frequency uj [§6.2.3 of 
iBinnev fc Tremaind l|1987f) ]: 

c.^ = «2-2.GE^+(^)^^ (3) 
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Fig. 1. — Toomre profiles Q{R) for various parameter sets. Panels (a) through (d) show the cflfcct of different choices for the spin 
parameter, A, and the disk mass fraction, (see J7|for a detailed description); each panel differs in only one parameter from those 
adjacent to it. Within each panel, profiles are shown for virial masses of 1 X 10* M0 (short-dashed line), 1 X lO'' Mq (dot-dashed), 1 X 10^° 
Mq (long-dashed), and 5 X 10^" Mq (solid). The Toomre parameter is temperature dependent: within each panel, dotted lines show where 
the Toomre threshold {Q = 1) would lie if the gas were at temperatures of lO** K (bottom), 300 K (middle), and 30 K (top). Disks are 
unstable below the Toomre threshold. To illustrate the relationship between Toomre stability and self-shielding, we also show the radius 
within which the gas is self-shielding (see JSJ for two different CBR spectral indices, rj = 0.7 (diamonds) and 2.0 (crosses). 



Here, k(R) is the epicyclic frequency of the orbit, and 
Cs(T) is the sound speed within the gas. Since the general 
waveform's time evolution is like e~*'^*, stability requires 
that cj^ > 0, and hence: 



(4) 



This is the Toomre stability criterion, written in terms 
of the Toomre parameter, Q. 

In order to determine the range of stable ^s, we com- 
bine equation ^ with the condition > to obtain: 



> 4 



(5) 



where the left-hand side is now recognisable as . Since 
> 0, stability is assured for: 



i > 4rit 



Also, for £^ -C ^crit' inequahty © reduces to: 



i < 



GS ' 



(6) 



(7) 



which is just the Jeans criterion for the gravothermal 
stability of a thin disk (see also Schaye 2004). These 
two cases can be interpreted as the zero pressure and 
zero rotation limits, respectively. Large perturbations 
are stabilised by the effects of rotation (inequality 0) , 



while small perturbations are stabilised by internal pres- 
sure (inequality |7| . 

Once the disk and halo have settled into dynamical 
equilibrium, all gross dynamical quantities (e.g. qh, S, 
y, k) become fixed. The only variable quantity in in- 
equality Q that can influence Q is then the sound speed, 
Cs — [kT / n)^/'^ . The amount of the disk that is unstable 
therefore depends strongly on T. This is illustrated in 
Figure m which shows profiles of Q{R) for the various 
parameter sets that we will consider, and the regions in 
which these disks would be unstable at temperatures of 
10* K, 300 K, and 30 K. 

As mentioned in ^ for low mass disks, thermal sta- 
bility at IC* K implies and is implied by gravothermal 
stability; conversely, low mass disks only become Toomre 
unstable at temperatures that are only accessible via the 
H2 cooling channel. Once H2 cooling is initiated, the gas 
makes a rapid tr ansition to a cold phase with T < 300 K 
ijNorman fc Soa ans 1993), which case it is virtually 
guaranteed (at least for M200 ^ 10^ M©) to be Toomre 
unstable. For these reasons, we use the H2 cooling rate of 
the gas — in chemical equilibrium and at T = 10"* K — in 
conjunction with the Toomre criterion for gravothermal 
instability to define our criterion for SF. 

Specifically, at a given point, star formation is deemed 
to occur wherever the H2 cooling rate exceeds the total 
rate of photoheating, provided that the gas be Toomre 
unstable at T > 300 K. Figure ^ illustrates the relation 



On Star Formation and Dark Galaxies 



5 



between these two criteria: gravothermally unstable re- 
gions lie below the Toomre threshold; we also show the 
radius within which the gas is self-shielding. 

4. THE MODEL — II 
THE DISTRIBUTION OF MATTER 

In the next three sections, we provide a detailed, tech- 
nical discussion of our three-step modelling procedure. 
The means by which we determine the distribution of 
matter in the protogalactic disk is described in this sec- 
tion. We discuss the solution for the equilibrium chem- 
ical composition of the gas in Sj^l and the calculation of 
the global minimum SFR is discussed in A reader 
wishing to avoid such a technical discussion may choose 
to skip to Sj3 

The solution for the volume distribution of the gas 
involves two separate calculations. First, the analytic 
model of Mo, Mao & White (1998; hereafter MMW) is 
used to determine the radial distribution of matter in the 
model protogalaxies. The MMW model has been shown 
to predict some global characteristics {e.g. size distribu- 
tion, velocity profiles, TuUy-Fisher relation) of nearby 
disk galaxies and damped Lyman-a systems, based on 
more or less generic predictions of A''-body cosmologi- 
cal simulations (MMW). Secondly, once the radial mass 
distribution is known, the vertical distribution of mass 
in the dis k is dete r mined usi ng an extension of t he for - 
malism of lBahcalll l|1984)) and lDove fc ThronsonI l)1993j) . 
Thes e two solution processes are described in H4.ll and 
M.'2\ respectively. 

In this section, and throughout the rest of this work, 
the symbols R and Z are used to represent the (cylindri- 
cal) galactocentric radius and the vertical distance from 
the midplane, respectively; the lowercase r is reserved for 
spherical coordinates. 

4.1. The MMW Solution for the 
Radial Mass Distribution of a Protogalactic Disk 

The MMW model operates in the 'adiabatic limit', in 
which the galax y formation process conserves an 'adia- 
batic invariant' l|Barnes fc White 1984). The final con- 
figuration is thus made independent of the assembly pro- 
cess, making it possible, for a known initial mass dis- 
tribution, to determine the final configuration, provided 
the form of the final disk mass distribution is specified. 

Following MMW, the dark and baryonic matter are as- 
sumed to follow the same initial mass distribution; specif- 
ically, the fitting formula proposed by Navarro, Frenk & 
White (1997; hereafter NFW): 



Ph{r) 



Pcrit 



(8) 



r/rh (1 -t- r/ni)2 

which has been shown to describe equilibrated halos in 
ACDM cosmological simulations. 

The halo scale length, rh, is defined in relation to the 
quantity r2oo, which is chosen such that the mean den- 
sity within a sphere centred on the halo and with ra- 
dius r2oo is 200 times the critical density for closure, 
Pcrit(^) = 3H\z)/8ttG; le.: 



Afaoo = M(r2oo) = 200 pcrit :j7r 



(9) 



Numerical experiments show r2oo to approximately de- 
marcate between regions in which material is still in- 



falling, and those in whi ch material tends tow ards equi- 
fibrated (radial) orbits l|Cole fc LacevllT99^ : r2oo and 
M200 are therefore referred to as the virial radius and 
mass, respectively. [SeefBullock e t al.l l|2001br) for a more 
detailed discussion of the relation between r2oo/-^^20o 
and the virial radius/mass.] The precise relation be- 
tween r2oo and r^ is set by the concentration parameter, 
c — r2QQli'h- The normalisation factor 5c in equation JSJ 
is defined as: 



200 3 , , , 
-c3^1n(l + c)-^ 



(10) 



Secondly, motivated b y the observed distribution of 
light in nearby galaxies ljFreemanlll970|) . the assembled 
gaseous disk is assumed to have an exponential surface 
density profile, S(i?): 



y(j}\ - ^^200 R/R^ 



(11) 



where the normalisation factor is chosen so that the total 
disk mass is a parametric fraction, nid, of the total dark- 
plus-baryonic mass. Thus, for a given M200 and rud, 
S(i?) is completely specified by a scale length, Rd. 

The iterative MMW solution for Rd relates the total 
energy of the pre-assembly halo, E, to the angular mo- 
mentum of the assembled disk, Jd, using the definition 
of the spin parameter. A: 



J' \E\ 
G2 Mi,, 



(12) 



where Jd — jdJ is a parametric fraction, jd, of the to- 
tal angular momentum of the pre-assembly, dark-plus- 
baryonic matter halo, J. The final result is the following 
expression for Rd (MMW): 



Rd = (2 fc) /fl r-200 
rud 



(13) 



where /c(c) and /i? (A, c, rud, jd) are dimensionless factors 
of order unity, dependent on the exact shape of the initial 
and final mass distributions, respectively. Wherever di- 
rectly comparable, our results are ind istinguishable from 
those of MMW and of lSchavel ll2fiol . 

4.2. The Vertical Mass Distribution 
of the Assembled Disk 

Once S(i?) is known, at a fixed R, the vertical den- 
sity profile for the disk gas, pd{R, Z) is governed by two 
equations — the Poisson equation: 



92$ 



AttG {Qh + Pd) 



(14) 



and the equation of hydrostatic equilibrium (in the Z 
direction) : 

2 dpd ^ 9$ 
-dZ^P'dZ 







(15) 



In these equations, $(i?, Z) is the gravitational poten- 
tial, ^r{R, Z) is the radial component of the potential, 
Qfi (r) is the density of the dark matter halo (determined 
numerically as in Appendix B), and a is the (thermal) 
gas velocity dispersion. Since the gas is assumed to be 
isothermal, a is taken to be constant throughout the disk. 
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Fig. 2. — Mass distribution cross-sections for M200 = 10® Mq, illustrating the importance of the quantity 5. Contours are of constant 
density, in units of mp cm~^. Each panel is SfJd on a side. The disk scale height Zd(R) is shown by the long-dashed line; profiles are 
computed out to 4Z;j. Models where g is computed as described in i|4.2l are shown on the right of each frame. For comparison, models with 
a lSoitzer. L1942D sech^ profile are shown on the left. 



Using the equation of hydrostatic equilibrium is equiv- 
alent to using the first moment in ITz of the coUisionless 
Boltzmann equation, with the velocity dispersion tensor, 
(see, e.g., §4 of Binney & Tremaine, 1987), set to 
5ija'^ . In the context of non-interacting stellar dynamics 
(c/. an equilibrated gas), this assumption would intro- 
duce an error of order {Z Zd/RRd) (Bahcall 1984). 

The disk is assumed to be in centrifugal balance in 
the R direction, and particles to travel stable, circular 
orbits; viz. V$fl = V^/R, where V{R) is the orbital 
velocity of the gas at a radius R. With this substitution 
in equation ifT^ . these two equations can be combined 
by differentiating equation (1 1511 with respect to Z, to give 



by dill 
((Bahc 



Bahcall 1984; Dove & Thronson 1993,) : 



Pd dZ 
AttG 



1 ^ 

Pd dZ^ 

AttG 



Pd 



1 



Qh 



AttG R OR ' ■ ^^^^ 

Following iBahcalll (|1984[) and iDove fc ThronsonI 

(^93), equation (|16|l is made dimensionless by defining 
a scale height for the disk Zd{R): 

( <r^ 

Zd{R)^[ ^^ , (17) 



2ttGpo{R)J 

where po{R) = Pd{R, 0) is the volume density of the 
disk gas on the plane, at a distance R from the galac- 
tic cen tre. [Note tha t Zd is twice the exponential scale 
height l)Bahcal]|ll98'^ .] This leads (again, for a fixed R) 



to an ordinary differential equation for the dimensionless 
density y{x) — pd{R, Z)/ po{R), where x — ZjZd'. 



2i(R,Z)y, (18) 



subject to the boundary conditions: 



y{x = 0) = 1 



and —— 



dx 



= 



(19) 



It is worth making a few comments here about the 
quantity ^ in equation (|18() . which is defined as: 



i(R,Z) = 



iirGpo 



AirGghir) 



1 dV^jR) 
R dR 



(20) 



but which can be rewritten in a more convenient form 
f ApDendix lB|) . If ^ is neglected altogether, then equation 
lfTH|) is analytic, yielding y{x) — sech^ x ljSpitzeil ll942). 
This formalism was first derived in the context of stel- 
lar dynamics close to the galactic plane and at modest 
distances from the galactic centre, in which case both V 
and Qh vary slowly across the volume of interest. Au- 
thors therefore typically neglect the dV^/dR term and 
treat Qh as constant {e.g. Maloney 1993; Dove & ShuU 
1994). 

This work, however, considers lower mass galaxies with 
manifestly non-flat rotation curves. Moreover, we are 
interested primarily in the outermost gas, since this is 
where the effects of the CBR are greatest; this is also 
where the effect of a non-zero ^ is most pronounced. It is 
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therefore inappropriate to adopt these standard approx- 
imations here. The effect of a non-zero ^ is ihustrated in 
Figure 121 which shows several model gas distributions as 
well as the equivalent distribution assuming y — sech^ x. 

In summary, the modelling procedure is as follows: 
y{x) is determined by numerical integration of equation 
(|18l) : po is then fixed with reference to the definition of 
viz.: 

Po = , . (21) 

a2 {2j-^dxy{R,x)f 

The presence of po in the definition of ^, as well as the 
implicit dependence of ^ on make it necessary to solve 
for y{x) iteratively: a trial value of po is used to deter- 
mine the value of ^ in the integration of equation (|18() : 
the value of the integral is then used in equation H21|l 
to obtain the next po- This process rarely requires more 
than several iterations for convergence. 

5. THE MODEL — III 
THE CHEMICAL COMPOSITION 
OF THE GAS AND THE CBR 

Once the distribution of the disk gas is known, we de- 
termine its equilibrium chemical structure. Nine chemi- 
cal species are identified within the gas, interacting via 21 
collisional and nine photodestructive processes, as out- 
lined in t|5.II The gas is subject to a UV — X-ray CBR 
field; we describe the CBR parameterisation, including 
gas self-shielding, in t l5.2l The solution for the equilib- 
rium chemical composition of the gas is discussed in ii5.3l 

5.1. Chemical Species and Processes 

Followin g Haiman, Thoul & Loeb (1996; hereafter 
HTL) and iAbel et ^ (fr99l . we distinguish nine chemi- 
cal species within the gas: H°, H^, H+, Hj, HJ, He°, 
He"*", He"'"'', and e~, assuming a (primordial) helium 
mass fraction of 0.24. These chemical species are al- 
lowed to interact via the 21 collisional processes identified 
by HTL. In addition, the nine photoionisation and pho- 
todissociation processes listed in Abel et al. ( 199^ are 
included. For ionisation o f He° and H° , the more recent 
cross- sections given by lYan. Sadeehpour fc Dalgarnol 
l|I998j) are employed. These expressions give much 
lower cross-sections at high energies than those used by 
^bcl ct al. (1997), which reduces the rate of heating in 
marginally shielded areas, but increases the penetration 
depth. The rate coefficients and cross-sections for all of 
these processes are listed in full in Appendix 1X1 

We completely neglect the effects of any dust or metals 
present in the gas, even though the IGM is expected to 
have seen conside rable enrichment l|Aguirre et al.ll200H 
iSchave el al.l2000|) . Metals can make a significant contri- 
bution to the cooling rate within the self-shielded region 



I|na,l1i fc Pallairmg^TKaFz. Weinberg fc HerTinuistjrmMh : 

moreover, both dust and metals act as sites for H2 pro- 
duction. These effects only become dominant, however, 
once two-body collisional production of H2 becomes in- 
efficient, well below the assumed temper ature of 10"* K 
ijBuch fc ZhanglI99HlGani fc Pallall99^ . Further, both 
dust and metals can absorb UV CBR radiation, so con- 
tributing to the level of self-shielding within the cloud. 
The neglect of dust and metals therefore minimises the 
predicted level of H2 production and cooling. 



For each individual chemical species, an expression can 
be written for the net rate of production/consumption 
in terms of the number densities of all nine species. 
It is most convenient and instructive to cast these ex- 
pressions in terms of the relative chemical abundances, 
Xi = rii/nn, where rii is the number density of species i, 
and rifj is the number density of hydrogen nuclei, viz.: 



E E 



dxi 
nH dt 

y y 

(22) 

Here, ap,q is the temperature dependent rate coefficient 
for the reaction between p and q, Tp is the rate of pho- 

toionisation/photodissociation of p, and bpq, dp^ — 0, 
±1, or ±2, depending on the number of i particles pro- 
duced/consumed in the reaction. Note that in the high 
density/zero radiation limit, the equilibrium abundances 
are independent of density; they are determined by the 
temperature dependent as alone. Equation H22|l encap- 
sulates a particularly 'stiff' set of nine highly nonlinear 
coupled differential equations for the chemical composi- 
tion of the gas at a fixed point (HTL) . 

5.2. The Photodestructive Effects of the CBR 

We consider the CBR in the range 
13.6 e\ <hv < 40 keV {ie. 912 A > A > 0.31 A— 
wc find that photoionisation of He and He+ by soft 
X-rays is important in H2 production, as well as being 
a significant photoheating mechanism.) Within this 
range, the CBR spectrum, Jcbk{^), is described using 
a simple power law, characterised by its spectral index, 
77, and a normalisation constant, J21, chosen such that: 



-^cbr(i', R, Z) = J21 



Hz 



^-t{u;R,Z) 



xlO"^' erg-' s~' Hz"' cm"^ Sr-V (23) 

Here, z^hi is the frequency of the Lyman-a photon: 
hvm ~ 13.6 eV. This parameterisation is the simplest 
possible; to assess its relevance, we direct the reader to 
the compendium of CBR observation s acros s the w hole 
electromagnetic spectrum provided by 'Henry' fl999f). We 
will trial two extreme values of 77 (see t i7.3(l , in an attempt 
to bracket the real situation. 
The frequency dependent optical depth, t, defined as: 



t{v- R,Z) = Y, 



Oiiiy) / dZ' ni{R,Z') 



(24) 



is used to modulate the CBR spectrum as it penetrates 
into the cloud, so mimicking radiative transfer. Here, 
aj(i/) is the frequency dependent cross-section for pho- 
todcstruction of species i, as given in Appendix 1X1 Note 
that this treatment of radiative transfer is 'one way' — 
no attempt has been made to account for radiation that 
crosses the midplane — leading to a very minor under- 
estimation of the photoheating rate close to the plane. 

We have made three substantial omissions in writing 
equations (|23|l and (|24|l . First, any extinction of the 
CBR due to the presence of dust or metals has been 
neglected. Secondly, secondary ionisations by energetic 
electrons from X-ray ionisations are ignored. While these 
ionisations have no effect on the heating rate, the en- 
hanced electron abundance would increase the H2 pro- 
duction rate (and the total rate of cooling); the omission 
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of secondary ionisations will therefore again lead to a 
lower limit on the shielded fraction. Thirdly, and most 
importantly, diffuse radiation produced by the gas itself 
has been ignored; that is, any photons (re-)radiated by 
the gas are assumed to be free to escape unhindered. 
This assumption is justified for the optically-thin outer- 
most gas; its validity within the self-shielded region is 
certainly problematic, but a detailed treatment of radia- 
tive diffusion through the cloud would be dependent on 
both the thermal and dynamic evolution of the disk dur- 
ing its formation. As has been stated, this is beyond the 
scope of the present undertaking. 

For each photoreaction except Solomon photodis- 
sociation of H2 (see below), the photodissocia- 
tion/photoionisation rate, F^, is then given by: 

F, - 27r / a,{v) , (25) 

where Vi is the threshold energy for dissociation of species 
i, and the factor of 2tt arises from integration over all 
angles of incidence onto the flat disk. 

Photodissociation of H2 occurs primarily in the 
Lyman- Werner b ands (11.26 — 13 .6 cV) via the two-step 
Solomon process (IStecher fc Wilhams 196'^ : 

H° +hv^ ^ 2H° . (26) 

This process is harder to deal with than other pho- 
toprocesses, since the structure of the molecular en- 
ergy levels (and so the cross-section for photodestruc- 
tion) is so much more complex. Moreover, Doppler 
broadening causes individual lines to overlap. Follow- 
ing ^^el ct al. (1997), absorption is therefore restricted 
to the Lyman band (12.24—13.51 eV); the CBR spec- 
trum is treated as flat in this narrow range. Thus, 
rather than explicitly considering H2 dissociating pho- 
tons longward of 912 A, the H2 photodissociation rate is 
determined with reference to the mean number of pho- 
to ns near 1000 A. The penetration pro bability formalism 
of lde Jong. Dalgarno fc Bolandl l)198f]t> . which modulates 
the unextinguished CBR field with a shielding factor, /3, 
is used to simulate H2 self-shielding, including the effects 
of Doppler broadening; the adopted expression for /3(th2 ) 
is given in Appendix 1X1 

5.3. Chemical Composition in Equilibrium 

Once the radiation density at a given point is known 
(and hence the rate coefficients for the various photo- 
chemical processes), it is possible to find the equilib- 
rium chemical composition at that point by solving the 
system of equations represented by equation H22|l . with 
the time derivatives on the left-hand side set to zero. 
We do this using a modified Netwon-Raphson algorithm 
(§9 of Press et al. 1992), where we compute the Ja- 
cobi matrix analytically. For dxi/dt = 0, only six of 
the nine equations remain independent; it is also nec- 
essary to ensure that the total number of hydrogen nu- 
clei, helium nuclei, and electrons are conserved separately 
l|Katz. Weinberg fc Hernnuistlll996l) . 

As a diagnostic for the solution procedure only, Figure 
Elshows the equilibrium chemical composition (in the ab- 
sence of photodestruction) as a function of temperature. 
In this plot, the lines show our steady-state solutions for 
the full set of equations encapsulated by equation (|22() : 



for comparison, the points show algebraic solutions for 
a reduced system of six equations obtained by omitting 
H~, H2, and H^. Note that these results differ signif- 
icantly from those in the directly comparable Figure 2 
of HTL, which were obtained by direct integration of 
equation H22|l until apparent convergence. In particu- 
lar, our predicted 11° abundance is one to two orders of 
magnitude lower than that of HTL, depending on the 
reaction rate expressions used (see Appendix El • The 
disagreement arises from the fact that the results shown 
in HTL had not yet fully converged (Zoltan Haiman, pri- 
vate communication — note that this does not affect the 
conclusions of HTL), as can be readily verified by direct 
substitution. 

It should also b e note d that, contrary to 
iHaiman. Rees fc Loehl l|1996D . assuming chemical 
equilibrium is expected to lead to a lower limit on the 
actual H2 abundance at T = 10'' K, and by extension 
on the H2 cooling rate. The assumed rate coefficient for 
H2 production is temperature independent, while that 
for consumption decreases with temperature; thus, the 
effect of a small decrease in temperature is an increase 
in the net ra t e of H-;. production [see, e.g., Figure 4 
of lAbel et"al] l)1997f) ]. The net rate of H2 production, 
therefore, can only go to zero if the H2 abundance is 
dropping (or if the gas temperature is rising). The H2 
abundance in an equilibrating gas will therefore always 
reach equilibrium from above. (For a time-evolved 
demonstration of this point, see Figure 4 of Haiman, 
Abel and Rees 2000.) 

Moreover, the situation of chemical equilibrium for 
H2 is physically improbable for temperatures between 
~ 10^ — lO'' K, because in the absence of a source of 
photoheating, the gas will be thermally unstable; gas is 
far more likely to cool past these temperatures without 
equilibrating. In this sense, we emphasise the limited 




Fig. 3. — Equilibrium abundances of nine chemical species in 
the absence of any photodestruction. Chemical species are as la- 
belled. The points show the equilibrium solutions for the reduced 
set of equations in which H~, Hj, and are omitted. All abun- 
dances are scaled relative to the total hydrogen nucleus number 
density, rajj- Note that the electron abundance closely tracks the 
proton (H"*") abundance. This plot is provided only as a diagnostic 
for our solution procedur e, an d should not be mistaken for actual 
evolutionary tracks, (see >l5.3i . 
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physical relevance of Figure 13 which is provided only as 
a diagnostic, and refer the reader to plots of nonequilib- 
rium abundances in cooling and collapsing gas clouds like 
those f ound in Figures 3— 5 of Oh & Haiman ( 200 2^), Fig- 
ure 3 o f lAbel et alJ l|1997i) , or Figure 6 o flShapiro fc Kand 
lfT987t> . 

6. THE MODEL ~ IV 
THE MINIMUM SFR 

At each point where H2 cooling is found to be suffi- 
cient to initiate gravothermal instability, we introduce 
an additional internal radiation field and determine the 
minimum ISRF required to preserve thermal balance at 
10^ K; it is then possible to relate this ISRF to a mini- 
mum global SFR, TZmin- In other words, we solve for the 
minimum ISRF-regulated SFR within the cloud. 

This is done by including a tenth equation in the sys- 
tem of equations described in ti5.3l requiring that the 
collisional H2 cooling rate (Ahs , given in Appendix 
be balanced by the net rate of photohcating, viz.: 

=0. (27) 

i 

Here, Q is the rate of photoheating due to absorption by 
species i, which (for all processes except Solomon disso- 
ciation of H2 — see Appendix A) is given by: 

t/(i^) 

df {hv - hvi) -J— ai{v) , (28) 

where the radiation flux J = Jqbr + cu^ now has two 
components. The first is due to the CBR, where Jcbr 
is defined by equation H23|l : the second is the new diffuse 
ISRF, characterised by the diffuse radiation energy den- 
sity, Ud- Note that we use this same J in equation 
so that the ISRF maintains thermal balance by some 
combination of photoheating and H2 photodissociation. 

The spectral shape of Ud is fixed with reference 
to the stellar popula tion modelling utility Starburst99 
ijLeitherer et al.iri999(l ^. which gives the SED for contin- 
uous star-formation, Lgg^u), with units of erg Hz^^ 
(1 Mq yr-i)-i. We adopt a standard SF scenario, which 
assumes a Salpeter IMF over a mass range of 1 — 100 M0, 
and a metallicity of 0.020 (Z©); our results are not 
strongly dependent on the choice of metallicity. Over 
the spectral range of interest, this SED converges after 
< 10 Myr (a few times the lifetime of O — B stars); we 
use this converged SED. It should be noted that the Star- 
burst99 SED only extends to 100 A; since self-shielding 
against the diffuse ISRF is ignored, this will lead to only 
a minor underestimate of the rate of Hel ionisation, and 
so the rate of heating. 

To convert this luminosity to a radiation density, it 
is necessary to make an estimate of the mean lifetime of 
each photon, i^, which is done with reference to the mean 
free path of a photon with energy hi/, i^, = {na,^)~^ = 
tvC. With these substitutions, Ud can be written: 

Udiv) = U Lgg{i^) ^c'Y^riiaiiv)^ . (29) 

^ Specifically, we used the data shown in Figure 2b, available at 
|http://www.stsci.edu/science/starburst997| 



While the spectral shape of Ud is given by L99, its inten- 
sity is thus characterised by a scale factor, U, which has 
the units of a SFR per unit volume. 

U is thus taken as the tenth unknown in a system of 
ten equations, representing the minimum SFR (per unit 
volume, and assuming chemical equilibrium) required to 
produce a sufficient ISRF to maintain thermal balance at 
T = 10'* K, and so restore gravothermal balance. Note 
that U gives the SFR required to ensure the stability 
of the volume dV, not the SFR within dV. Figure 0] 
shows the resultant ISRF field in cross-section for several 
exemplary models. A global SFR, T^min, is then found 
by integrating U = dTZmin/dV over the whole disk. 

Clearly, this is not a solution for the exact problem. 
That said, T^-min is expected to provide a reasonable esti- 
mate for the minimum SFR in the cloud, since 'maximum 
damage' has been allowed; photons are delivered directly 
to wherever they are required. Under these conditions, 
so long as the actual SFR exceeds T^-min, the non-star- 
forming remainder of the disk will remain stable; other- 
wise, some part of the disk will become unstable and the 
SFR will increase as a result. 

7. THE MODEL — V 
CHOICE OF PARAMETERS 

The model now outlined consists of three components 
described by nine parameters. In this section, the fiducial 
values adopted for these parameters are discussed. 

Cosmology enters our calculations via Pait, which is 
used to determine r2oo for a given M200 (see equa- 
tion ini; we adopt the now st andard WMAP cosmol- 
ogy (Ben ne tt et al. 2003; Socrg el et all 120031. Further, 
model disks are assumed to form at zero redshift. Since 
Pcrit decreases with decreasing redshift, r2oo, and hence 
Rd, will be greater for smaller disk formation redshifts 
(see equation 9); at a fixed R/ Rd, this minimises pd and 
so self-shielding and H2 production, in keeping with the 
ethos of the model. 

7.1. The Progenitor Halo: M200; c, A 

The maximum virial mass considered is M200 = 5 x 
10^* Mq, typical of a nearby LSB galaxy. The minimum 
mass, M200 ~ 10^ Mq, is set by the form of the ana- 
lytic model: below this limit, the solution process for Zd 
and po fails to converge at R ^ Rd- Mathematically, 
the long, slow rise of the rotation curve for lower mass 
galaxies makes ^ a large negative number at these radii, 
which eventually drives the vertical density profile, y{x), 
to increase with distance off the plane (see equation EOJ. 
In this unphysical case, we truncate the integration, and 
treat the stationary point as a sharp boundary. It may 
be that lower mass halos do not form exponential disks. 

The concentration parameter, c, is chosen using 
a simple scaling relation dr awn from Figure 6 of 
iNavarro. Frenk fc Whitd l(T99^ : 

where AI^ is the (cosmology dependent) 'nonlinear' mass, 
as defined in that work. 

Numerical experiments show that the spin distribution 
in hierarchical clustering scenarios is well described by 
a 'log-normal' distribution centred at Aq — 0.04, and 
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A = 0.04; ma = 0.05 
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Fig. 4. — Contours of Ud, shown in cross section. Ud is the minimum ISRF energy density (at 1000 A) required for thermal balance, 
determined as per Contours are shown (in units of the Habing flux) for example model (A, m^) = (0.04, 0.05); clockwise from top-left 
are models with M200 = 5 X 10^", 10^", 10®, and 10* Mo- Note the different ranges of the vertical axes in the upper {Rd) and lower {iRd) 
panels. We only solve for Ud at radii where the gas is Toomre unstable at T > 300 K (see Figure 0. Overlaid are contours of njj (dotted 
lines), in units of cm""^. 



with w idth 5\ = 0.5 (jWarren et al.lll992l: iBullock et all 
l2001aj) . Two values for the spin parameter, A, are 
adopted to sample this distribution: a 'typical' value of 
0.04, and a 'high' value of 0.10. These values approxi- 
mately correspond to the 50 and 90 percent points of the 
expected distribution (sec also MMW). 

7.2. The Disk: ma, jd> c, T 

Ba sed solely on the WMAP cosmology ijSpergel et al.1 
I2003j) . baryons are expected to make up 16 percent (by 
mass) of all matter. However, not all of the baryonic ma- 
terial initially within an overdense region will ultimately 
settle into the disk. Moreover, of the mass that does 
collect at the centre of the potential well, some fraction 
may be lost, e.g. thro ugh evaporation after reionisation 
ijShaviv fc Deke]|l200^ . Two fiducial values for the disk 
mass fraction are therefore adopted, nid — 0.05 and 0.10 
(see also MMW). 

Assuming that the halo's angular momentum origi- 
nates from tidal torquing, it is natural to assume that 
both the dark and baryonic components obtain the same 
specific angular momentum (Fall & Efstathiou 198(|); 
thus, at the time of virialisation, = j^. While 
A^-body and hydrodynamic simulations predict that 

i d/rud is significantlv less than unity for assembled d isks, 
Navarro fc Whitell994t lN avarro fc Steinmetzll997D . the 
MMW model requires that jd/rud be close to unity in 
order to produce disk sizes and properties that are in 
agreement with observations. Accordingly, we assume 
jd = TTid for the assembled disk. 



The velocity dispersion, ct, is assumed to be constant, 
with the value of 8.2 km s~^; this corresponds to the 
thermal velocity dispersion at T = lO'' K. Our chosen 
value is larger than the canonical value of 6 km s^^, 
which, using a recalibrated Toomre criterion (Q < 2), 
has been shown to predict the location of star forma- 
tion thresholds in around 2/3 of nearby spi ral galaxies 
l|Kennicutljll989t iMartin fc Kennicuttll200lD . However, 
iSchavc ( 2004) has recently suggested that this agreement 
is largely coincidental: for typical HSB disk galaxies only, 
the point where the disk is Toomre unstable assuming 
Cs = a = 6 km s~i fortuitously coincides with where the 
transition to the cold phase becomes possible. 

7.3. The CBR Field: J21, ry 

The CBR parameters are fixed using the value of 
the low-redshif t Hi photoionisation rate, Fhi, obtained 
bv iScott et"aH (2002 ). They give a best fit value of 
Thi = 2.0^11 X 10"" s-i for z < 1, where Tm is related 
to the UV — X-ray CBR spectrum by: 



Thi = 47r 



dv — — — aHi(j^) . (31) 



Here, ohi is the known cross-section of Hi (see Appendix 
For a power law spectrum, this integral can be per- 
formed analytically, which makes it possible to fix the 
normalisation constant J21 for a given spectral index 77. 

We consider two fiducial values for 77. The first, rj = 2, 
mimics a typical quasar spectrum (|Scott et al. 2002): the 
second, 77 — 0.7 is chosen to match cosmological simula- 
tions of the Lyman-a forest and the observed X-ray CBR 
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up to 40 keV llHaiman. Rees fc LoebllTOQ^ . Using the 
iScott et aJ] l)20fl2tl value for Fm. these two choices corre- 
spond to values for J21 of 0.091 and 0.071, respectively. 

8. RESULTS 

We determine the amount and distribution of gas (as- 
suming chemo-dynamical equilibrium at T = 10*^ K) that 
is subject to H2 cooling in the presence of the CBR. 
Where the cooling rate due to H2 processes exceeds the 
total rate of heating due to the CBR (ze. where the gas is 
thermally unstable at 10"* K) the gas is deemed 'shielded'. 
If this gas is also Toomre unstable at T > 300 K {ie. if it 
will be gravothermally unstable in the cold phase), then 
it is deemed 'unstable'; this gas is then prone to star for- 
mation. Additionally, by determining the ISRF intensity 
required to restore thermal balance at T = 10"* K at each 
'unstable' point (so restoring the situation of gravother- 
mal stability), we estimate the minimum global SFR for 
each model. 

A few trial models have been recalculated with greater 
precision to ensure that numerical effects do not signif- 
icantly alter our results; our principal results, shown in 
Figures and [7| are convergent to better than 1 percent. 

8.1. Localised Properties of Model Disks 

In Figure |5l we present some of the 'local' scaling re- 
lations that arise within the models. Within each panel, 
the two distinct curves in regions of marginal shielding 
correspond to the two choices of CBR spectral index: the 
upper curve to rj = 2.0, the lower one to 77 = 0.7. 

First, in panel (a), the mass fraction of shielded gas is 
plotted as a function of surface density. These results are 
in good agreement with the first estimate made in ^2.'2\ 
A universal vertical density profile would ensure that all 
columns with the same S would be self-similar; the scat- 
ter in this plot — particularly in the rj = 2.0 case, which 
has a higher proportional X-ray flux — therefore demon- 
strates of the importance of detailed modelling of the 
vertical gas distribution. Within the scatter, the trend 
with increasing mass is right to left; more massive disks 
are better able to self-shield. 

In panel (b), the minimum ISRF intensity required to 
maintain thermal balance at 10^ K, Ud, is plotted as a 
function of gas density. The ISRF is presented in units of 
the Habing 'flux' (the observed local ISRF energy density 

A-i cm-3 = 1.33 X 



(a) Self-Shielding 



at 1000 A), = 40 X 10"^*^ erg 



10~29 erg Hz~^ cm"^ ( Habing ' 1968H . Within the self- 
shielded region, Ud = 0.16 n^^^ g; this relation has been 
overplotted. 

Finally, in panel (c), the minimum SFR (per unit 
projected area) required to preserve thermal balance, 
T^min/ dA, is plotted as a function of Hi column den- 
sity. Again, within the self-shielded region, the curves 
converge towards a power law: TZnun/ dA « 1.1 x 
10-" (TVm/lO^i cm-2)3 = 2.2 x 10"" S^j. The re- 
su lts plotted in this pan el agre e very well with those 
of iCorbelli. Galli fc Fallal l|1997D . Note that surround- 
ing the completely self-shielded region is one in which 
Tlmin/ dA exceeds this power law: in this region, the CBR 
actually aids SF, as the increased electron abundance due 
to X-ray ionisation of Hel enhances H2 production. 

Since thermal balance within the self-shielded region 
is achieved primarily through photodissociation of H2 
(rather than heating by the CBR), the results of the 
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Fig. 5. — Local scaling relations arising in the model: (a) The 
fraction (by mass) of gas that is shielded from the CBR as a func- 
tion of surface density; (b) the minimum ISRF intensity required 
to maintain thermal balance at 10'* K as a function of density; 
(c) the minimum SFR required to maintain thermal balance as a 
function of HI column density. Representative points are drawn 
from each of the 32 parameter sets we have modelled. Within each 
panel, the two distinct curves correspond to r; = 2.0 (upper curves) 
and r] = 0.7 (lower curves). Panel (c) should not be interpreted in 
terms of a Schmidt relation. The SFR plotted in panel (c) is the 
SFR required elsewhere in the cloud to prevent SF in a column of 
density A^hi; it is not the SFR expected in that column. 
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Fig. 6. — The fractional mass that is unstable to SF in the absence of an internal ISRF. The left panel shows models with r\ = 0.7; 
the right shows r] = 2.0. Within each panel, from top to bottom, parameter sets are (A, m^) = (0.04, 0.10), (0.04, 0.05), (0.10, 0.10), 
(0.10, 0.05). In other words, lower spin galaxies are more unstable than higher spin; then for the same spin, a higher fractional disk mass 
reduces stability. Note that for an exponential disk, only ~ 90 percent of the mass resides within AR^i; this point is shown by the dotted 
line. Even for the 'worst case' CBR, the majority of baryonic matter is unstable to SF in at least 50 percent of all galaxies with baryonic 
masses greater than a few times 10® M©; the same is true for all galaxies with baryonic masses greater than ~ 10^ Mq. 



ISRF calculation are sensitive to the H2 photodissocia- 
tive flux alone. The predictions for Ud are therefore ro- 
bust, in that they should not be sensitive to the assumed 
ISRF spectral shape, or to the exact distribution of the 
gas. However, the link between and TZ^m/dV is more 
tenuous, since the exact relation depends on the extent, 
location, and character of SF within the cloud. 

8.2. Integrated Properties of Model Disks 

Figure El shows the fractional mass of gas available for 
star formation for each of the 32 parameter sets that we 
have modelled. All else being equal, high spin galaxies 
are less able to self-shield than those with lower spin. In 
the high spin case, Rd is typically 8 — 10 times greater 
than for low spin, which reduces the surface density at a 
fixed R/Rd by a factor of order 50 — 100. 

Also, a CBR spectrum with a shallower spectral in- 
dex has a greater proportional X-ray flux, and so is able 
to penetrate deeper into the cloud. At a fixed point, 
this increases the rate of photoheating; at the same time, 
the free electrons from helium ionisations lead to an in- 
creased H2 production rate. Thus, in the r] = 0.7 case, 
even though the CBR aids SF near the shielding bound- 
ary, this boundary is pushed deeper within the gas. 

Recalling that the adopted values for A represent the 
50 and 90 percent points of the expected distribution, 
the majority of the disk is predicted to be unstable to 
SF in at least 50 percent of putative dark galaxies with 
baryonic masses greater than 5 x 10^ M©. For baryonic 
masses greater than 10^ M©, the majority of mass is 
unstable in as much as 90 percent of all protogalaxies, 
depending on the distribution of and the true CBR 
spectrum. 

This result is directly contrary to that of 



IVerde. Oh fc Jimenej lj2002f) . who have predicted 
that all disks in halos with M200 < 10^ M© should 
remain dark. However, these authors have assumed 
a constant Cs = cr = 6 km s~^ when calculating the 
Toomre parameter, which implies T k, 8500 K; in the 
context of this model, this is tantamount to ignoring H2 
cooling. 

The predicted global minimum SFR within the pu- 
tative dark galaxies, 7?.min, is plotted in Figure [T] as a 
function of the Hi mass. A/hi; in ah cases, A/Hi/™dAf200 
is between 0.56 and 0.70. For comparison, some ob- 
served values of Mhi and SFR for actual dwarf galax- 
ies have been overplotted. We have included: a sam- 
ple of six "relatively nearby" {D < 15 Mpc) Im and 
seven "more distant" (30 <D < 100 Mpc ) Im a nd Sm 
galaxies from iHunter. Elmegreen fc Baker! l)1998(l ; four 
Hl-bright nearby galaxies, two each with "low" (< 10~* 
M© yr-^ a nd "high" (> 10"^ M© yr^i) SFRs, from 
lYoung et all (|2003); se ven LSB d warf galaxies and four 
"normal" dwarfs from Ivan Zee et al. (1997) ; 121 irreg- 
ular g alaxies from the sample of Hunter fc ElmegreenI 
lj2004(l . which span more than 8 mag in absolute magni- 
tude and surface brightness; and ESO215-G7009, an ex- 
tremely Hl-rich dwarf galaxy, with the highest accurately 
measured Hi massdight ratio to d ate, Mm/ Lb ~ 20 
ijWarren. Jerien fc Koribalski 2004*). 

For each combination of A and nid, T^min oc M^^. We 
note that this relation has the same power as the em- 
pirical Schmidt SF law, SFR/dA cx Ej^f (Schmidt 1959; 
Kennicutt 1989, 1998; but see also Boissier et al. 2003). 
Moreover, even though these samples were selected to in- 
clude widely disparate populations, these populations are 
not readily distinguished in Figure [3 for Mhi ^10® M©, 
nor is the exceptional ESO215-G7009 clearly differenti- 
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Fig. 7. — Minimum SFR required for stability, 7?.min, plotted against HI mass, Mm. As in Figure|H] like models are connected with a 
dashed hue to guide the eye: from top to bottom, models are (A, m^) = (0.04, 0.10), (0.04, 0.05), (0.10, 0.10), (0.10, 0.05); t] = 0.7 (2.0) 
models are plotted with crosses (diamonds). Note that the top and bottom pairs of lines represent the 50 and 90 percent points of the 
expected spin distribution, respectively, and that we have deliberately attempted to place a firm lower bound on T^-min- For comparison, 
observed val ues of Mht a.nd SFR for a variety of dwarf galaxies are overplotted, including: 'nearby' (large pluses) and 'distant' (small pluses) 
dwarfs from Hun ter. Elmegreen fc Bakerl 119981k 'low' SFR (small crosses) and 'high' SFR (large crosses) nearby Hl-bright dwarfs from 
lYoung ct al. |2003); 'LSB' (small diamonds) and 'normal' (large diamonds) dwar fs from [v an Zee ct al. (1997); 121 members of a sample 
of irregular galaxies (points), spanni ng a large range in galactic parame ters, f rom [Hunter &c Elmegreen Ii2004) : and the extremely Hl-rich 
dwarf ESO215-G7009 (large square) IWarren. Jerien KoribalskM2004l) . See il9.4l for a discussion of the significance of this relation. 



ated from the others. 

These results are essentially independent of the CBR 
spectral index. For the Hi fraction, which is determined 
principally by the outermost gas in the cloud, this is 
because the CBR spectrum has been normalised using 
Phi- On the other hand, 7<^min is determined principally 
by the dense gas closest to the plane of the disk and to the 
centre of the galaxy, where the effect of the CBR is least. 
In combination with the results in FigureEl this gives the 
surprising result that T^min is essentially independent of 
the amount of gas that is deemed 'unstable'; it depends 
on Mhi alone. 

9. DISCUSSION 
9.1. The Effects of Substructure and Turbulence 

It is difficult to be quantitative about the effect of sub- 
structure, but, qualitatively, self-shielding is increased 
where the surface density is enhanced. Since, to first 
order, the CBR prevents a skin of fixed column depth 
from entering the cold phase, removing matter from 
a marginally shielded column and adding it to a well 
shielded column will increase the overall shielded frac- 
tion. That is, the omission of substructure such as bars 
and turbulent cloudlets should again lead to an underes- 
tim ate of the amount of gas pron e to SF. 

iWada. Meurer fc NormanI l)2002^) have shown that sta- 
tistically stable turbulence, leading to a multi-phase in- 



terstellar medium, can be produced and sustained with 
no stellar energy input. It is conceivable that this tur- 
bulence could dominate the velocity dispersion, so pro- 
viding a third mechanism for support against gravita- 
tional collapse. The turbulent velocity dispersion re- 
quired to support gas at T 300 K, however, is typi- 
cally more than 15 km s~^, whereas in the model these 
authors present, the turbulent velocity dispersion is only 

1 — 5 km s~^.^ 

Moreover, iWada fc NormanI l|200lD have shown that 
the timescale for cooling is far shorter than the several 
dynamical times required to establish this turbulent mo- 
tion. In their model, the central gas cools within 0.1 Myr, 
whereas turbulent effects first become apparent only af- 
ter ~ 2 Myr, and stabilise after ~ 10 Myr. While tur- 
bulent effects are not unexpected in these galaxies, we 
therefore expect them to be insufficient to stabilise the 
gas, and to emerge only after the first bout of star for- 
mation. 

9.2. An Important Sanity Check 

The crux of our argument is the assumption that, if 
and when H2 cooling is initiated, the gas will rapidly 
evolve into a cold phase at T ~ 300 K. In W2.'6\ we have 
argued that this transition will be fast in comparison to 
the dynamical timescale; we are now in a position to 
verify that this transition will, in fact, lead to Toomre 
instability. 
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Fig. 8. — Minimum temperature achieved by cooling gas, Tmin, 
plotted against the critical temperature below which the gas be- 
comes Toomre unstable, Tq, for representative points in regions 
deemed unstable. Where Tq > Tniin, gas is able to make an evo- 
lutionary transit to a Toomre unstable cold phase, in line with the 
argument presented in JS] This is essentially always true. 



Detailed integration of the coupled rate equations en- 
capsulated by equation H22() have found that gas cooling 
from T > 10** K to T < 300 K dev elops a non-eq uilibrium 
abundance of Hy , xho "•^ 10~'^ llShapiro fc Kanslll98^ 
ISusa et al.ill998t lOh fc Haim'all2002|) . nearlv indeoen- 
dently of the initial conditions. This occurs at T 
4000 K, when (two-body) coUisional H2 production and 
consumption processes become inefficient, leaving an 
un-replenishing 'freeze-out' abundance of H2. In the 
presence of a photodestructive radiation field, these 
molecul es can be quickly des troyed, preventing further 
cooling. lOh fc HaimanI lj2002() have shown that, because 
of the scalings in the problem, the final temperature that 
a parcel of gas (at fixed n) is able to cool to depends only 
o n the ratio J/n . 

lOh fc Haiman (2002) describe a method for estimating 
the minimum temperature achieved by a cooling par- 
cel of gas, Tmin, by finding the temperature for which 
the timescales for cooling, tc, becomes longer than that 
for photodissociation, — ^n]- With the approxima- 
tion fx T^, and assuming the freeze-out abundance, 
xb_2 — 10^'^, this leads to the following expression for 

T„,i„ « 9000 (^^^) K . (32) 

(Recall that (3 is the H2 self-shielding factor, discussed 
in ^b.2\ ) It should be noted that this analysis does not 
include dust catalysed H2 production at T < 3000 K. 

In Figure IHl we plot T^in (for gas in unstable regions) 
as a function of the critical temperature below which the 
gas becomes Toomre unstable, Tq. Points that lie above 
the line Ty^in = Tq represent regions where the gas will 
transit to a Toomre stable cold phase, in violation of the 
assumption on which we base our SF criterion. There 
do appear to be very rare violations of this assumption. 



however these points are those closest to the shielding 
boundary (far from the midplane), in the most Toomre 
stable regions (far from the centre) and make an insignifi- 
cant contribution to the important integrated quantities. 

9.3. Implications for Star Formation 

Our model argues for a subtle but significant alteration 
to the (often unstated) conventional wisdom concerning 
SF: that Hi is the raw material from which stars form. 
Within this paradigm. Hi is commonly referred to as 'un- 
processed' gas, and it is common and reasonable to in- 
terpret the Schmidt 'law' as a causal relation between an 
amount of gas and the level of SF. 

The primary role of H2 cooling in this picture lies in 
the fragmentation of a collapsing gas cloud into proto- 
stellar cloudlets. We have shown, however, that H2 can 
also be important in initiating the first Toomre/ Jeans 
instability. This is particularly important in the case of 
low mass galaxies which would remain Toomre stable if 
their thermal evolution were governed by Hi alone. 

The importance of this slight revision lies in the inter- 
pretation of Hi observations. The picture outlined above 
interprets Hi as, in some sense, the 'original' state of hy- 
drogen. In the alternate scenario, however, H2 can be 
seen as a 'default' state for hydrogen: wherever possi- 
ble. Hi will always develop into H2. Moreover, since H2 
cooling virtually guarantees Toomre instability, some ex- 
ternal source of UV radiation is required to prevent the 
production of stars from within the gas by dissociating 
H2 into Hi. In this new paradigm. Hi is thus a product of 
SF rather than being its primary fuel, and the Schmidt 
relation is symptomatic of the mode of SF. 

In other words, the Schmidt relation should not be 
interpreted as 'given more Hi, more stars will form', but 
instead as 'more UV radiation (produced by short-lived 
O — B stars) is required to support a bigger Hi cl oud.' 
[For a detailed exposition of this argument, see lAllenI 
ll200l .] 

This revision is not so controversial: the impor- 
tance of H2 cooling in low mass galaxy formation has 
been well established by works such as Peebles fc Dicka 
^1968!), "LeiDD fc Shull (1984), Haiman. Thoul & Loe3 
((199(11 ■ and .Kepner. Babul fc Spereel (19971 Although 
these works are typically (but not exclusively) focused on 
the effect that radiative feedback from the first genera- 
tion of stars has on galaxy formation, they clearly demon- 
strate the importance of H2 in the global thermal and dy- 
namic evolution of galaxies before protostellar fragmen- 
tation. Closer to home, t he established physics o f pho- 
todissociation regions [see IHollenbach fc TielenFJ dTofll 
for a comprehensive review] demonstrate the importance 
of the UV ISRF in producing Hi in our galaxy. 

9.4. The Mbi—SFR Relation 
and DDO 154~like Galaxies 

Our model predicts a hitherto undocumented relation 
between Hi mass and SFR, which seems to hold over the 
range Mm ~ 10*— lO^" Mq: SFR oc M^^. Within the 
mo del, this relation is a manifestation of self- regulating 
SF ()Lin fc Murravll992tl0^ukai fc NishilT999|) . in which 
the ISM is kept warm and stable by a UV ISRF that is 
maintained by constant regeneration of O — B stars. 

The slope of the relation is driven by the distribution of 
the densest gas in galaxy centres, since it is this gas that 



On Star Formation and Dark Galaxies 



15 



requires the most radiation to maintain thermal stability. 
The scatter is set primarily by the distribution of galactic 
spins, as well as variations in the relative mass of the disk 
in comparison to the halo, since, for a fixed mass, these 
determine the density near the centre. 

Because the model was designed to place a lower limit 
on the SFR, it cannot predict the zero-point of such a 
relation. Within the self-shielded region, the H2 cooling 
rate is typically less than 1 — 10 percent of the net cool- 
ing rate; a significantly higher ISRF might be required 
to balance these other cooling mechanisms. Moreover, 
if there were, for example, a systematic run with mass 
in the relative contribution of metal-line cooling in com- 
parison to H2 cooling, the inclusion of metals might alter 
the predicted slope of the relation. We stress that this 
would only be true, however, if metal-line cooling were 
important in initiating the transition from warm to cold. 

The general agreement between the character of the 
predicted and observed Mhi — SFR relations lends weight 
to the physical picture on which this simple model is 
based: the results in Figure [7| suggest that the well ob- 
served population of dwarf galaxies represent the mini- 
mum rate of ISRF-regulated SF in galaxies. 

Moreover, in the range Mhi ~ 10* — 10^° Mq, it seems 
that the majority of dwarfs are forming stars in the same 
way: 'extreme' objects like DDO 154 and ESO215-G7009 
are not distinguished from 'normal' dwarfs in the Mhi — 
SFR plane. In response to the question "Why was so 
little gas processed into stars?", therefore, these galax- 
ies' high Hi mass:light ratios are not due to inefhciencies 
in the present mode of SF. Instead, they may be, for 
example, a reflection of late formation times or differing 
environmental effects. 

9.5. The Complete Prevention of SF 

Figures and [7| show that SF is completely prevented 
in only one of the models: (A/200, A, nid) — (10* M0, 
0.10, 0.05). As is apparent in Figure ^ this is due 
primarily to rotational support rather than an inabil- 
ity to self-shield: the gas is self-shielding at i? « 1 (2.5) 
Rd for 7] = 0.7 (2.0), but is never Toomre unstable at 
T = 300 K. However, if the gas were able to cool to 
~ 30 K {ie. if it could initiate metal line cooling), then 
up to 70 percent of the gas would become available for 
SF, and we would predict TZmin ^ 10~^ M0 yr^^. 

9.6. The Detectability of Dark Galaxies 

Throughout this work, since the focus has been on 
gravothermal stability in the absence of internal SF, it 
has been implicitly assumed that the gas is predomi- 
nantly Hi (and not Hll or H2). Clearly, a very large stel- 
lar component would be necessary to completely ionise 
the gas, and it has been shown that the gas becomes 
gravothermally unstable in the presence of small amounts 
of H2. It has also been argued that, once SF begins, the 
ISRF will regulate the SFR by dissociating H2 into Hi. 

Hi is detectable by the hyperfine transition at ~ 21 cm, 
provided that the excitation temperature, Tgpin, exceeds 
the mean background radiation temperature, Tbg ~ 3K. 
Tspin can be written (Kulkarni & Hcilcs 1988): 



where, neglecting Hi excitation by electron collisions, 
y ~ (T/1000 K) / (nHi/0.2 cm'^). From Figure it 
is apparent that nni is essentially always in the range 
10~^ — 10^ cm~^ within the self-shielded region, in which 
case Tspin = 800 — 9800 K. Moreover, since Hi is so opti- 
cally thick to Lyman-a, even the low predicted SFRs act 
to thermalise Tgpin via scattered Lyman-a photons. 

Without exception, all models are thus predicted to 
have detectable Hi emission, independently of their op- 
tical properties. Joint optical-plus--21cm surveys are 
therefore predicted to yield complete samples of galaxies. 

10. THE EXISTENCE OF DARK GALAXIES — 
SUMMARY 

We have developed a model for the long-time configu- 
ration of a baryonic dark galaxy, in order to determine 
whether such a galaxy, having formed, could conceivably 
remain dark. The model has the following features: 

• The gas is in a dynamically equilibrated disk (ro- 
tationally supported in the radial direction and in 
hydrostatic equilibrium vertically) at the end of the 
galaxy formation process. 

• The gas is isothermal at T = 10^ K, which repre- 
sents the endpoint of its thermal evolution in the 
absence of H2 cooling. Together, the assumptions 
listed above limit the model to considering galaxies 
with M200 > 10* Mq. 

• The gas is in chemical equilibrium, in the presence 
of a UV— X-ray GBR. The GBR can prevent SF 
by preventing the gas from undergoing a transition 
(driven by H2 cooling) from a warm, largely stable 
phase to a cold, mostly unstable phase. 

We evaluate the stability of this situation of ther- 
mal, chemical and dynamical equilibrium by determin- 
ing whether or not H2 cooling can induce gravothermal 
instability, and hence SF, in the face of the GBR. Ad- 
ditionally, wherever an ISRF is required (in addition to 
the GBR) to counterbalance H2 cooling, we place an ap- 
proximate lower bound on the expected SFR. 

We find that, in the absence of an internal UV radia- 
tion field, gas that is sufficiently self-shielded to support 
an H2 abundance as low as 10""* is unable to remain 
gravothermally stable. This statement remains true inde- 
pendently of any other cooling processes, including metal 
line cooling: so as to place a firm lower bound on the 
amount of gas that will become gravothermally unsta- 
ble, we have ignored the significant shielding and cooling 
effects of dust and metals when calculating the H2 cool- 
ing rate at 10"* K. Moreover, if the gas is even mildly 
enriched, then its final temperature d rops by an order of 
magnitude from ^ 300 K to ~ 30 K ijNorman fc Spaan^ 
Fr997); particularly for the models most affected by the 
GBR, this could significantly reduce stability against SF 
(Figure HI). 

Our conclusions can be summarised as follows: 

1. The model predicts a correlation between Hi mass 
and SFR, which is observed in a wide variety of 
dwarf galaxies (see Figure [71). This suggests that 
dwarf galaxies represent the minimum levels of 
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ISRF-regulated SF in the universe. Further the- 
oretical and observational studies of this relation 
should prove fruitful. 

2. With only one exception, the modelled galaxies 
cannot avoid SF indefinitely, and so will develop 
a stellar component. Even for baryonic masses as 
low 5 X 10^ Mq, the majority of gas in greater than 
50 percent of all halos is unstable to SF in the ab- 
sence of an internal ISRF (Figure EJ. 

3. Without exception, all models would have de- 
tectable Hi emission. That is, whatever dark or 
dim galaxies exist should be detectable in HiPASS- 
class surveys. 



galactic Hi clouds in the Local Group that did not also 
have a stellar component. Sim ilar results are available 
for ot her groups and clusters (jZwaanll200lUWaugh et alJ 
120021) . Within these limits, there is no reason to suspect 
a population of undetected baryonic dark galaxies. 
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Above the detection limit of Mri ~ 7 x 10*^ Mq 
l|Mever et all 12004,1 . HiPASS did not detect any extra- 
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APPENDIX 

A. CHEMICAL PROCESSES IN THE DISK 

The full list of reactions considered in this study is given in Table Al, along with their rate coefficients. All quantities 
in this section are listed in cgs units. The symbol r„ is used to denote T~". Note that the rate coefficient for reaction 
(21) has two terms; the second term arising from dielectric recombination ljCenl ll992'l. The rates given for reactions 
(16) and (17) both are valid in the low density {nn < cm~'^) l imit. 

Also note that three of the expressions given in Table 3 A of .Haiman. Thoul fc Loebl l)1996ri differ from those given 
in the works that those authors cite. They give the exponent of T3 in t he expression for the ra t e coef ficient of reaction 
(2) as 0.2; it should read -0.2 [This typographical error is not present in lHaiman. Rees fc Loebl l|1996|) ]. More seriously, 
the multiplying factor in the expression for reaction rat e (3), given as 5.57, shoti l d read 3.23. Also, the exponent of T 
in reaction rate (6), which should be 2.17, is omitted in lHaiman. Thoul fc Loebl l)1996j) . These alterations reduce the 
equilibrium abundances of H2 by several orders. 



# 




Reaction 


Rate Coefficient, a (cm'^ s ^) 


Ref. 


(1) 


Hn 


' e~ — 


* H+4 2e- 


5.85 X 10-11 rl/2 e-157809.1/T (1+7.1/^)-! 


Cen92 


(2) 


H+4 


' e~ — 


♦ H"4 hu 


8.40 X 10-11 r-1/2 T^°-^ (1 + r0-7)-i 


Cen92 


(3) 


HO4 


' e~ — 


> H-4 hu 


3.23 X 10-1^ Ti/2 


HTL96 


(4) 


Hg4 


- e~ — 


* H-4 HO 


2.7 X 10-* g-43000/T 


HTL96 


(5) 


HO4 


H-- 


* HO4 e- 


1.3 X 10-^ 


Raw88 


(6) 


HO4 


H-- 


* 2HO4 e" 


5.3 X 10-20 T2-17 e-8750/T 


S&K87 


(7) 


e- + 


H-- 


* HO4 2e- 


4 X 10-12 T e-8750/T 


S&K87 


(8) 


H+4 


H-- 


* 2H0 


7 X 10-'' T-1/2 


D&L87 


(9) 


H+4 


H-- 


. H+4 e- 


4 X 104 T-i--! e-15100/T : T > lO* K 
10-* . r < lo* K 


S&K87 


(10) 


H+4 


RO- 


* H+4 hu 


3.255 X 10-20 T - 4.152 x 10-25 - 6.157e-l' : T > 4000 K 
1.38 X 10-23 2-1.845 . 200 < T < 4000 K 


RDB93 


(11) 


H+4 


H|- 

h|- 


. H+4 H" 


2.4 X 10-0 e-21200/T 


S&K87 


(12) 


HO4 


. HO4 H+ 


6.4 X 10-10 


HTL96 


(13) 


H-4 


> HO4 H" 


5 X 10-*^ T-1/2 


D&L87 


(14) 


e-4 




. 2H0 


1.68 X 10-* (r/300)-0-29 


HTL96 


(15) 


HO4 HO- 


* 3H0 


6.11 X 10-14 e-29300/T . J, > 7390 ^ 












2.67 X 10-15 e-(6''50/T)2 ^ ^ ^ ^^gp ^ 


L&S83 


(16) 


HO4 HO- 


^ H!^4 2H0 


5.22 X 10-14 e32200/T . y > 7291 K 










3.17 X 10-15 e-(4060/T)-(7500/T)2 . y < ^^91 K 


S&K87 


(17) 


HO4 


- e~ — 


^ 2HO4 e- 


4 38 X 10-10 TO'35 e-102000/T 


S&K87 


(18) 


He04 


- e~ — 


* He+4 2e- 


2.38 X 10-11 Tl/2 e-285335.4/T (1+7.1/2)-! 


HTL96 


(19) 


He+4 


- e~ — 


^ He++4 2e- 


5.68 X 10-12 Tl/2 g-631515.0/T (i + T^l/S^-l 


Cen92 


(20) 


He+4 


- e~ — 


^ He'' 4 hu 


1 50 X 10~"^^ 2^—0.6353 










4 1.9 X 10-3 T-3/2 e-4™000/T(i 3g-94000/T) 


Cen92 


(21) He++4 


- e~ — 


* He+4 hu 


3.36 X 10-10 T-1/2 T3-0-2(i + TgO-')-! 


Cen92 



Table Al_— -Reaction rate c oefficients. The e xpressions given for reactions (16) and ( 1 7) ar e valid for rifi < 10^ cm J^. References: 
(Ccn92)'Ccn' (1992); (D&L87) Dalgarno & Lcdd (1987); (HTL96) Haiman, Thoul & Locb (1996); ('LfcS83') ITepp fc ShulJI (19831: (Raw88) 
itewlings (1988); (RDB93) Rawlings. Drew & Barlow (1993); (S&K87) Shapiro & Rang (1987). All coefficients are given in cgs units. Note 
that the rate coefficients for reactions (2), (3) and (6) given in lHaiman. Thoul fe Loebl 11996f) disagree with those given in the works those 
authors cite (see text). 



For photoionisation and photodissociation reactions, we employ the reaction cross-sections given by lAbel et all 
()1997|) . except for ionisation of He" and H", for which we use the more recent expressions given by 
lYan. Sadeghpour fc Dalgarnd (^98). The different asymptotic behaviour of the expression for He" ionisation has 
a significant effect on the chemical structure of our models, but does not drastically affect our conclusions. When 
calculating the cross-section for H2 ionisation, we assume all H2 to be in the para- configuration. This is appropriate 
only when reactions (5) and (12) are the dominant H2 production mechanisms; within the model these reaction rates 
are typically seven orders of magnitude higher than for reaction (13). 

To account for self-shielding of H2 against LW photons, including the effects of Doppler broadening and line overlap, 
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we use the shielding factor of Ide Jong. Dalgarno fc Bolandl l)1980() : 

1/2- 

/ Tn. \ 

X erfc 




-1/2 / - ^ 1/2N 



where erf c denotes the co mplimentary error function, and d = 9.2 x 10"'^ a^^ cm km s""'^ is called the Voigt parameter. 
Following Ide Jonj l)1977|) . we assume that each S olomon dissociation heats the gas by 2.5eV. 
Finally, we use the lengthy expression given bv lLepp fc ShulM l)1983j) for the rate of H2 coohng, Ahs: 

LvL = 8.18 X 10"^'^ (^CKiQ riHO + an rifjo 
LrH = 3.90 X 10-19 e-6"8/^ 

LrL = (l.38 X 10-22 ^-9243/T^ 

Q„ = 1.2 (nHo)°■^^ + (nHof^^ 

where each of the Ls has units of ergs s~^. As given here, this expression is valid only for t emperatures abo ve 4031 
K. For the calculation presented in ^9.21 we use: ~ 4 x 10""^^ T^, which is a fit given bv lOmukai fc Nishil fl9991 
and is a good appr oximation for 600 5, T < 3000 K and n < 10^ cm-'^. For all other cooling processes, we use the 
expressions given in lHaiman. Thoul fc Loebri)1996tl . and references therein. 

B. CALCULATING THE QUANTITY ^ 

The quantity ^(i?, Z) in equation (|16|l is related to what iBahcall l|1984>) calls the effective halo potential, except 
that any vertical variation in ^ is neglected in that work. The effect of explicitly including 1^ in determining the volume 
density of the disk gas is illustrated in Figure [3 In this section, we motivate and describe the method by which we 
have calculated ^ in our integration of equation H16() . 

Throughout this section, we use the symbols R and Z to represent the (cylindrical) galactocentric radius and the 
vertical distance from the midplane, respectively; we reserve the lowercase r for spherical coordinates. We will also 
use the symbols M, p, and V here to refer to quantities which can be expressed analytically, as distinct from M, g, 
and V, which must be computed analytically. 

As described in M4.2I ^ is defined as: 

Here, {R, Z) is a sum in quadrature of contributions due to the (spherical) dark matter halo and the (cylindrical) 
baryonic disk, viz. ljFreemanlll970t iBinnev fc Tremaindll987j) : 

V^R, Z) = Vlir) + V!{R) 



GMh{r) , GmrfA/200 , 2(r ,..^r. . ,..^r. ,..^\ . R 



• 2y'^h{y)K^{y)-h{v)K^{y)j ■■ V ^ , (B2) 

where /„ and Kn are nth order modified Bessel functions of the first and second kinds, respectively. Note that there 
is an implicit thin disk assumption embedded in this expression, since it neglects any variation of Vd with Z. 

The final halo mass distributio n, Aih(r). is determine d using the method described by MMW, which assumes that 
the 'adiabatic invariant', rM{r) ijBarnes fc Whitdll984j) . is conserved throughout the galaxy formation process. That 
is, a test particle initially located at a mean radius moves through successive radii r in such a way that, at any given 
time: 

n M,{ri) ^ r M{r) ^ r [Mh{r) + M^R)] . (B3) 

Here Mi{ri) is the (assumed) initial dark-plus-baryonic mass profile, Md{R) is the (trial) final disk mass profile, and 
Mh{f) is the (unknown) adiabatically contracted halo profile, viz.: 

Mh{r) = {I - ma) M,{n) . (B4) 

Combining equations ljB3|l and ljB4|l . we obtain an expression which we solve numerically for as a function of r: 

Mrf(r) + (1 - mrf - ^) M,{n) = Q. (B5) 

The resultant can then substituted back into equation (|B4p and thence l|B2|) to determine V. 

Of course, determining dV^ /dr is doubly expensive, since this numerical procedure must be repeated twice in order 
to estimate the derivative; the same is true of Qh'- 

1 dMh{r) 

^hKr) = — J — ^ • (B6) 

47rr^ ar 
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However, we can rewrite the halo contribution to the dV^/dr term in equation I|B1|I as ijKuiiken &: Gilmorelll989() : 



dR 



^ 9r _9 \ GMh{r) 
R OR [ r 
]_ R [G dMhjr) 
R r \ r dr 
1 fG 



GMh{r) 



in r Qh{r) - (1 - rn^) 



GAhin) 



AttG Qhir) - (1 -rud) 



GAUin 



where we have used equations (|B6|) and (|lj4() to get from the s econd to the third line. 
When this expression is substituted back into equation ljBl|) . the Qh terms cancel, and we are left with: 



iirGpo 



(1 - md) 



GMhin) 



dR 



(B7) 



(B8) 



In this new expression, both Mh and dV^/dR (jBinnev fc TremaindlToSTjl are known analytically. 

Thus, by incorporating a lookup table for ri{r), we have developed a means of retaining all contributing terms in 
equation l|l()|l . without significantly increasing computation time. Moreover, with our rewritten definition of ^, we 
have eliminated the need to compute dAih/dr and dV^ /dr numerically. We foimd these modifications to decrease the 
computation time for the vertical distribution of the gas by a factor of roughly one hundred. 



